Mechanical control of spin-orbit splitting in GaAs and InGaAs epilayers 
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Time- resolved Kerr rotation spectroscopy as a function of pump-probe distance, voltage and 
magnetic field is used to measure the momentum-dependent spin splitting energies in GaAs and 
InGaAs epilayers. The strain of the samples can be reproducibly controlled in the cryostat using 
three- and four-point bending applied with a mechanical vise. We find that the magnitude of the 
spin splitting increases linearly with applied tension and voltage. A strain-drift-diffusion model is 
used to relate the magnitude of the measured spin-orbit splitting to the amount of strain in the 
sample. 

PACS numbers: 71.70.Ej, 71.70.Fk, 72.25.Dc, 72.25.Rb 



Potential applications in spintronics^ and quantum 
information processing^ rely upon an understanding of 
the effect of electric fields and strain on electron spins. 
Strain reduces the symmetry of a crystal, which intro- 
duces momentum k-linear terms to the Dresselhaus^ and 
Bychkov-Rashba^ spin splittings. These strain-induced 
effective magnetic fields can be used to generate elec- 
tron spin polarization electrically^ and coherently ma- 
nipulate spins using electric fields and in the absence of 
magnetic fields^, but they also contribute to more effi- 
cient spin relaxation^. In addition, recent steady-state 
measurements^'^ have shown that the spatial period of 
strain-induced spin precession is independent of the ap- 
plied electric field, which demonstrates the robustness of 
strain-induced spin precession for applications in func- 
tional spin-based devices. 

Here we employ mechanical three- and four-point 
bending to tune the tensile strain of GaAs and In- 
GaAs epilayers while performing low-temperature time- 
resolved magneto-optical spectroscopy to determine the 
magnitude of the strain-induced spin splitting. The 
samples are contacted so that an in-plane electric field 
can be applied to impart an average drift velocity 
to the optically-excited electron spins. Kerr rota- 
tion measurements as a function of magnetic field and 
pump-probe distance are performed for different ap- 
plied electric fields, and we observe that the spin split- 
ting increases with increasing drift velocity and tensile 
strain. Unlike previous measurements that introduced 
strain through heterostructure engineering and lattice- 
mismatched growth^, these measurements are able to 
map out the strain dependence in a single sample and 
without the complications of strain relaxation. The vise 
geometry allows for repeatable tensioning of samples and 
precise control over the strain level. 

The samples are grown using molecular beam epitaxy 
on semi-insulating (001) GaAs substrates. We examine 
both n-doped GaAs and InGaAs epilayers. The GaAs 
samples are comprised of 100 nm undoped GaAs buffer 
layer, 400 nm Alo.7Gao.3As, and a 500 nm Si-doped GaAs 
epilayer. Samples with carrier densities of 2 x 10^^ cm~^ 
and 4 x 10^^ cm~^ were measured, but since they exhibit 



qualitatively similar behavior, we show only data for the 
2 X 10^^ cm~^ n-GaAs sample in this paper. The InGaAs 
sample is composed of 300 nm of growth-interrupted 
GaAs buffer layer, 500 nm of Si-doped Ino.04Gao.96As 
with a carrier concentration of 3 x 10^^ cm~^, and 100 
nm of undoped GaAs. The lattice constant of the InGaAs 
layer is matched to that of GaAs, as confirmed using x- 
ray diffraction. Since InGaAs has a larger natural lattice 
constant than GaAs, the InGaAs layer is compressively 
strained in-plane as grown. 

The samples are patterned into mesas using pho- 
tolithography and a chemical wet etch and then con- 
tacted using annealed Ni/AuGe [Fig. 1(a)]. The channels 
have a width = 120 /im and length / = 310 /im between 
the contacts and are aligned such that an electric field E 
= V// can be applied along either the [110] {x) or [110] 
{y) directions. The samples are then mounted into either 
a three- or four-point mechanical bending vise. While 
the four-point bending geometry produces strain that is 
uniform between the central two point s^^, higher maxi- 
mum values of strain can be achieved in the three-point 
geometry before structural failure occurs. In the case of 
the four-point bending, the mesas are patterned far from 
the four contact points, minimizing local strain varia- 
tions. The sample is cooled to a temperature T = 60 
K, at which nuclear polarization is negligible^^. For the 
field-dependent measurements, a magnetic field is applied 
along X. A cross-section of the measurement geometry is 
shown in Fig. 1(b). The vise is tightened along the z 
direction, introducing tensile strain along x. The optical 
measurements are performed along z and measure the 
spin polarization along [001]. 

Time-resolved Kerr rotation spectroscopy^^ is used to 
monitor the electron spin dynamics in the samples. In 
this technique, a Ti: Sapphire laser, tuned to the ab- 
sorption edge of the material that we wish to probe 
(wavelength A = 818 nm for GaAs and A = 850 nm for 
Ino.04Gao.96As), produces a train of ~250 fs pulses at a 
repetition rate of 76 MHz, which are split into a pump 
beam (2 mW) and a probe beam (200 /iW) with 30 jiia 
diameters. The pump pulse is circularly-polarized and 
excites a spin-polarized electron population in the epi- 
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FIG. 1: (a) Patterned sample. Mesa is shown in grey, and 
metal contacts are shown in gold, (b) Measurement geometry 
showing orientation of sample, mechanical vise, electric and 
magnetic fields and optical measurement, (c) (top) Kerr rota- 
tion at temperature T = 60 K on 2 x 10^^ cm~^ GaAs epilayer 
as a function of appUed magnetic field for electric field E = 32 
V cm~^ and pump-probe separation d = 38 /xm at time delay 
At = 13 ns for strained states 2 (black), 3 (red), 4 (green) 
and 5 (blue), (bottom) Kerr rotation on InGaAs epilayer for 
unstrained (black) and strained (red) states. 




-40 -20 20 40 
B{mT) 



-40 -20 20 40 
B{mT) 



FIG. 2: Kerr rotation at temperature T = 60 K and time 
delay At = 13 ns as a function of applied magnetic field Bext 
and pump-probe separation d on 2 x 10^^ cm~^ GaAs epilayer 
sample for electric field E = 32 V cm~^ for (a) strained states 
2 (b) 3 (c) 4 and (d) 5. The color scale is the same for all four 
plots. 



layer. A linearly-polarized probe pulse is incident on the 
sample at time At later, which is controlled using a me- 
chanical delay line. The electron spin polarization in the 
sample is measured by detecting the change in the po- 
larization axis, or Kerr rotation, of the reflected probe 
beam. Time-resolved measurements at temperature T = 
60 K and an applied magnetic field B = 0.2 T show that 
the electron g-factor is -0.43 for the GaAs epilayer and - 
0.48 for the InGaAs epilayer and that the transverse spin 
coherence time is 40 ns for GaAs and 3 ns for the InGaAs 
sample. For the spatially-resolved measurements, a step- 
per motor-driven mirror changes the spatial separation 
d between the pump and probe beams The applied 
electric field causes the spins to drift with an average ve- 
locity Vd and imparts a non-zero average momentum (k) 
to the electron spin packet. 

In order to determine the k-dependent internal mag- 
netic field Bint, we measure Kerr rotation as a function 
of the applied magnetic field Bext at At = 13 ns for var- 
ious E and pump-probe distances. The presence of Bint 
modifies the symmetric cosinusoidal field-dependent sig- 
n3]kiM^. When Bint is along the same direction as Bext, 
the signal becomes centered about -Bint, but if Bint and 
Bext are perpendicular, as is the case when the electric 
field and Bext are applied along the same direction, the 
data can be fit to the equation: 

KR = Acos{gf,B^/m^~+Bl',At/h) (1) 

where A is the amplitude, g the effective g-factor of the 
sample, jub the Bohr magneton, and h is Planck's con- 
stant over 27r. For E = 32 V cm~^ and At = 13 ns, 
the center of the electron spin packet is observed to be 



at pump-probe separation d = 38 /im. Figure 1(c) (top) 
shows data for the GaAs sample for different amounts 
of tension applied using the four-point bending vise. We 
measure Kerr rotation in the channel for the unstrained 
case and for five increasing levels of strain, which we la- 
bel as strained states 1-5. The vise is tightened by 
the same amount between each of these strained states. 
Increasing the strain in the GaAs sample decreases the 
signal amplitude and increases the spin precession fre- 
quency and Bint- The change in the amplitude is due to 
a decrease in the spin lifetime, which we confirm using 
time-resolved Kerr rotation. In contrast, measurements 
on InGaAs, shown in Fig. 1(c) (bottom), show that the 
application of tensile strain via the three-point bending 
vise leads to some relaxation of the compressive strain 
due to lattice mismatch. The Kerr rotation amplitude 
is increased, and Bint is decreased. Measurements were 
also performed on the 2 x 10^^ cm~^ GaAs sample for 
channels where E was applied along y and perpendicular 
to Bext and the direction of the tensile strain; in this ge- 
ometry. Bint is parallel to Bext, and the values obtained 
for Bint were 22% smaller. 

We explore the strain-dependent spin-orbit splitting in 
the GaAs sample as a function of electric field and pump- 
probe distance. In Figures 2(a)-(d), we show Kerr rota- 
tion as a function of applied magnetic field and pump- 
probe separation for strained states 2, 3, 4, and 5 for E 
= 32 V cm~^ and At = 13 ns. The color scale is the 
same for all four plots. Again, we observe that the spin 
precession period and amplitude decrease with strain and 
that an increase in Bint lowers the central peak. In these 
measurements, we also observe the effect of spin diffu- 
sion, which is manifest in the spatial dependence of Bint 
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FIG. 3: (a) Effective magnetic field Bint and (b) amplitude 
A as a function of pump-probe separation d on the 2 x 10^^ 
cm~^ GaAs epilayer sample in strained state 4 for electric 
field E = 8 (black), 16 (red), 24 (green), 32 (blue), 40 (cyan) 
and 48 V cm~^ (magenta). Symbols are data, and fines are 
fits as described in text. 

and is due to the spread in drift velocities of the spin 
packet. 

In order to characterize the voltage-dependence of Bint, 
we fit the amplitude A as a function of d with a Gaus- 
sian function to determine the center position of the spin 
packet dc- We then use the results of a linear fit of Bint 
to determine the value of Bint at dc for each voltage. In 
Figure 3, we show the data (symbols) and fits (lines) for 
A and Bint foi" the GaAs sample for strained state 4 as a 
function of d and for various values of E. From these fits, 
we obtain the spin-splitting energy Aq = g /j^b Bint at 
the center of the spin packet as a function of Vd^ which is 
plotted in Fig. 4 for increasing amounts of strain in the 
2 X 10^^ cm~^ GaAs sample. As observed previously^, 
the data can be fit to a line, where the slope P = Ao/vd 
can be used to characterize the observed effect. We plot 
f3 for each of the strained states in the inset of Fig. 4 
and observe that /3 increases for increasing amounts of 
applied tension. In comparison, previous measurements^ 
on GaAs strained by the removal of the underlying sub- 
strate showed that P = 99 neV ns /im~^. Similar mea- 
surements of the InGaAs sample reveal that /3 = 72 neV 
ns /im~^ when unstrained and = 40 neV ns /im~^ when 
strained with the three-point bending vise. 

Although we tightened the mechanical vise by the same 
amount between each of the strained states, slip and play 
in the vise make it difficult to determine the amount of 
bending and strain in the sample by the mechanical dis- 
placement Az. In order to estimate the amount of strain 
for each of the strained states, we solve a strain-drift- 
diffusion model^'^^ to determine the spatial spin preces- 
sion period as a function of £. Although this model was 
developed for steady-state measurements, as described in 
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FIG. 4: Spin-splitting energy Ao as a function of drift ve- 
locity Vd for 2 X 10^^ cm~^ GaAs epilayer sample when un- 
strained (black symbols) and strained states 1 (red), 2 (blue), 
3 (green), 4 (orange), and 5 (purple). Linear fits are also 
shown. (inset) (3 [neV ns /im~^] for unstrained (0) and 
strained states 1-5. The dashed line indicates /3 from mea- 
surements on a GaAs membrane in Ref. 

Ref. 1^, using the spatial spin precession period^^ SPP 
= 2 TT h/ /3, we can relate P and e. The model parame- 
ters used are similar to those used in Ref. pfl5j with the 
exception of the spin diffusion constant D. A value of 
D = 283 V/cm^ was obtained from the measurements in 
Fig. 3. The strain-drift-diffusion equations are 

Oip^c = O2PZ (2) 

Oipy = Ospz (3) 

O^Pz + + Ospy = -G, (4) 
where the operators O1-O4 are defined as follows, 

Oi = D\/^ + /iE • V - {CsDf - - (5) 
O2 = -Cb, - Cs{2D^ + pE,) (6) 



Os = CB^-Cs{2D—^pEy) (7) 
O4 = DV^ + /iE • V - 2{CsDf - - (8) 

Ts 

Here pi is the {x,y,z} component of electron spin polar- 
ization, Gz is a Gaussian source function with a FWHM 
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FIG. 5: (a) Spin polarization pz as a function of distance 
solved using the strain-drift-diffusion model for e = (black), 
1.6 (red) and 2.2 x 10~^ (green), (b) The linear relation 
between f3 and strain e from the strain-drift-diffusion model 
(line). The symbols show (3 obtained from measurements on 
the 2 X 10^^ cm~^ GaAs epilayer sample, (c) Estimates of 
strain e for the strained states from the measured (3. 

of 30 /im, D is the spin diffusion constant, /i is the mobil- 
ity, Ts is the spin relaxation time, E is the applied electric 



field, B is the applied magnetic field, Cs = C3 m* e/h'^, 
hCsi = g M B^, m* is the electron effective mass, and 
e is the strain. We assume a value for the spin-strain 
coupling coefficient C3 = 4.0 eV A, as in Ref. 15]. The 
equations are solved over a 1 x 1 mm field using a finite 
element method. For the strain calibration, the external 
applied magnetic field was set to zero, and an electric 
field of 33 V cm~^ was used. We obtain solutions of the 
spin polarization as a function of position for varying val- 
ues of e between and 0.001. Figure 5(a) shows three 
line cuts of pz taken along the direction of the electric 
field for different solutions with e = 0, 1.6, and 2.2 x 
10~^. We fit these line cuts to determine the spatial spin 
precession frequency as a function of e. A linear fit of /3 
as a function of e yields a slope of 22.27 ± 0.96 in units 
of P per 10~^ unit strain. Figure 5(b) shows /3 as a func- 
tion of s. This relation allows us to assign strain values 
to all six tension states using their measured values of /3 
[Fig. 5(c)]. 

In summary, we have performed quantitative measure- 
ments of the spin-splitting energy as a function of volt- 
age and strain on samples mounted in a mechanical vise. 
Using a strain-drift-diffusion model, we are able to esti- 
mate the strain in the devices and calibrate the observed 
spin-splitting. This spin-splitting can be used to locally 
and coherently manipulate electron spins and electrically 
drive spin resonance^' 
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